Epitaxial growth with pulsed deposition: 
Submonolayer scaling and Villain instability 
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It has been observed experimentally that under certain conditions pulsed laser deposition (PLD) 
produces smoother surfaces than ordinary molecular beam epitaxy (MBE). So far the mechanism 
leading to the improved quality of surfaces in PLD is not yet fully understood. In the present work 
we investigate the physical properties of a simple model for PLD in which the transient mobility 
of adatoms and diffusion along edges is neglected. Analyzing the crossover from MBE to PLD, 
the scaling properties of the time-dependent nucleation density as well as the influence of Ehrlich- 
Schwoebel barriers we find that there is indeed a range of parameters where the surface quality in 
PLD is better than in MBE. However, since the improvement is weak and occurs only in a small 
range of parameters we conclude that deposition in pulses alone cannot explain the experimentally 
observed smoothness of PLD-grown surfaces. 
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I. INTRODUCTION 

Pulsed Laser Depostion (PLD) is an increasingly used 
growth method which is particularly suited for the fabri- 
cation of multilayer thin films |Q] . The basic mechanism 
of this technique is that the target material is ablated by 
a pulsed laser and then deposited in pulses on a substrate 
so that many atoms arrive at the surface simultaneously. 
Experimentally each pulse lasts for about a few nanosec- 
onds, and the time between two pulses is of the order of 
seconds. PLD differs significantly from ordinary molec- 
ular beam epitaxy (MBE), where atoms are deposited 
continuously. Growth experiments show that the sur- 
faces in both cases have different morphologies. In par- 
ticular it was observed that in some cases PLD leads to 
smoother surfaces than MBE However, as there are 
only few theoretical studies so far, the mechanism lead- 
ing to smoother surfaces in PLD is not yet understood 
so that it is impossible to predict the growth quality for 
various materials and growth parameters. Therefore the 
aim of the present work is to improve our understand- 
ing of PLD by investigating the growth morphology for 
different growth parameters in a simple model. 

Without going into experimental details it should be 
mentioned that the physical conditions for PLD are far 
less well defined than for MBE. The particles deposited 
may be atoms, clusters or even droplets. They may arrive 
with energies ranging from 0.1 eV to 1000 eV. If very en- 
ergetic particles arrive at the surface their kinetic energy 
is converted into heat at the surface, changing locally the 
effective mobility of the particles for a short time. In the 
present work this transient mobility is neglected, i.e., it 
is assumed that the particles arrive at the substrate sur- 
face with energies of about 0.1 eV which are of the same 
order of magnitude as in the case of thermal deposition. 

Recently Narhe et al. investigated the island statistics 
in the coalescence regime for a PLD process where tin 
droplets are deposited on a sapphire substrate H . It was 



observed that the scaling differs significantly from MBE 
due to the large fraction of multiple droplet coalescence 
under pulsed vapor delivery. While these results are valid 
for high deposition energies, the present work investigates 
a different physical regime, namely pulsed deposition of 
individual atoms at low energies and the formation of 
islands in two dimensions. 

The model investigated in the present work is defined 
as a solid-on-solid growth process on a square lattice of 
L X L sites with integer heights representing the configu- 
ration of the adsorbed layer. The particles are deposited 
in pulses with an intensity / which is defined as the num- 
ber of particles per unit area deposited per pulse. The 
duration of a pulse is assumed to be zero and the tran- 
sient enhancement of the mobility of freshly deposited 
adatoms is neglected. The model is controlled by three 
parameters, namely, the intensity / of the pulses, the 
diffusion constant D, and the average flux density of in- 
coming particles F. One of these parameters can be fixed 
by choosing the time scale so that we can use / and D/F 
as independent parameters. The dynamic rules are de- 
fined as follows, (i) In each pulse IL^ atoms are instan- 
taneously deposited at random positions on the surface, 
(ii) Between two pulses a time interval At — I/F elapses, 
in which adatoms diffuse on the surface with rate D. If 
Ehrlich-Schwoebel barriers are present the rate at which 
particles hop down the edge of an island is reduced, (iii) 
If two atoms at the same height occupy neighboring sites 
they stick together irreversibly, forming the nucleus of a 
new island or attaching to an already existing island. In 
our model the influence of edge diffusion is neglected so 
that islands grow in a fractal manner before they coa- 
lesce. We present results from kinetic Monte Carlo sim- 
ulations on a lattice of 400 x 400 sites. Finite size effects 
are discussed in the Appendix. 

The paper is organized as follows. In the following sec- 
tion we first recall the properties of PLD without Ehrlich- 
Schwoebel barriers and its crossover to MBE in the limit 
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MBE PLD 

1: Molecular beam epitaxy (left) compared to pulsed 



laser deposition (right) for D/F = W and / = 0.01. The 
figure shows typical configurations after deposition of 0.05 
ML. 



of very low intensities ||]. In Sec. [II we investigate the 
scaling behavior of the time-dependent nucleation den- 
sity in PLD, extending our previous analysis in Ref. 



In Sec. [V we study the influence of Ehrlich-Schwoebel 
barriers on the morphology of surfaces grown by PLD, 
finding a parameter range where layerwise growth is im- 
proved compared to MBE. 
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FIG. 2: Average island distance £d versus pulse intensity I 
for different values of D/F. The data points between the two 
dotted vertical lines are used to determine the exponent v. 
The dashed line is the corresponding power-law fit with the 
slope —v = —0.26(1). The inset shows the saturation levels 
of £d for small / as a function of D/F. The dashed line has 
slope 7 = 0.17(1). 



II. FROM MBE TO PLD: 

CROSSOVER AND SCALING 

If the intensity / is very small PLD behaves essentially 
in the same way as MBE. In fact, at the lowest possible 
intensity, where only one atom per pulse is deposited, 
MBE and PLD are equivalent up to minor statistical dif- 
ferences resulting from finite size effects (see Appendix). 
Using lattices of 400^ sites we therefore restrict our anal- 
ysis to the range D/F < 10^ where these finite size effects 
are negligible. 

At high intensities the growth morphology of PLD dif- 
fers significantly from MBE. As shown in Fig. |l], there 
are much more nucleations at an early stage, although 
the effective flux of incoming particles is the same in both 
cases. The two regimes are separated by a crossover at 
a certain intensity Jc, where the number of deposited 
atoms per pulse is of the same order of magnitude as the 
average adatom density in the corresponding MBE pro- 
cess. Obviously, if the pulse intensity is much higher than 
the MBE adatom density, the adatoms nucleate much 
faster forming many small islands. For this reason PLD 
is expected to yield more homogeneous surfaces, which 
explains the technological interest in this method. 

Let us now study the crossover from MBE to PLD in 
more detail ||] . The adatom density in MBE averaged 
over space and time is known to scale as {D / F)'^'^^^ [Q, 
where 



_ 1 _ J 1/6 for compact islands 
~ df+A ^ \ 0.18 for DLA fractals § 



Thus the crossover takes place at the critical intensity 

I, oc {D/F)^'<-\ (2) 

A quantity which distinguishes the two different growth 
modes shown in Fig. |l| is the average island distance. 
Performing numerical simulations we realized that the 
scaling regime of this quantity is not only restricted by 
finite-size effects but also by lattice effects, which be- 
come relevant for high intensities as well as for low val- 
ues of D/F, where the average distance of islands is of 
the order of a few lattice constants. Combining these 
bounds we find that for a system with 400^ sites finite- 
size and lattice effects are negligible in the parameter 
range 10^ < D/F < 10^ and / < 10"^ 

Fig. H shows the average island distance as a function of 
the intensity for various values of 13/ -F at a fixed coverage 
of 0.2 monolayers (ML). For low intensities I < Ic the 
island distance depends only on D/F. Plotting these 
saturation values versus D/F (shown in the inset of the 
figure) one recovers the well-known power law for the 
island distance in MBE 

Id oc {D/Fy (3) 

with 7 = 0.17(1). This estimate lies between the val- 
ues for compact growth and diffusion-limited aggrega- 
tion (see Eq. (0)), supporting the assumption that the 
islands are characterized by an effective fractal dimension 
1.6 <df <2. 

For intensities I > Ic the island distance is indepen- 
dent oi D/F and can be described by a power law 

£d cx r-" (4) 
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FIG. 3: Data collapse of the curves in Fig. | according to (|) FIG. 4: The nucleation density at the bottom layer versus 
using the exponents 7 = 0.17 and v = 0.26. The knee of the monolayer time during the deposition of one monolayer. The 
curves marks the crossover from MBE- to PLD-like behavior. dashed line has the slope 1. 



with an exponent v = 0.26(1). The independence of the 
island distance of DjF is a result of the high density of 
adatoms which nucleate so quickly that they do not make 
use of their full mobility given by D/F. 

Combining Eqs. (||) and (Q) we obtain the scaling 
form H] 

Id cc {D/Fr ■ Kl/h) , (5) 

where his a, scaling function with the asymptotic behav- 
ior 

As shown in Fig. this scaling form leads to a convincing 
data collapse. Moreover, the independence of D/F at 
high intensities together with Eq. (^) implies that 

— ^ ~ / ^-^^ compact islands 
^ ~ 1 - 27 ~ \ 0.28 for DLA fractals ' ^ ' 

confirming our numerical estimate v = 0.26(1). 



III. TIME-DEPENDENT NUCLEATION 
DENSITY 

In a recent paper we investigated the time- 
dependent nucleation density in PLD, reporting an un- 
usual type of scaling behavior. In order to avoid the 
crossover to MBE, we considered the limit — > 00 
where Ic — 0. In this limit all adatoms nucleate or attach 
to existing islands before the next pulse arrives. Fig. ^ 
shows a log-log plot of the the nucleation density at the 
bottom layer n{I , <d) as a function of the coverage Q — Ft 
for various intensities /. Obviously this quantity does not 
display ordinary power law scaling since it is impossible 
to collapse the curves by shifting them horizontally and 



vertically. However, in Ref. we observed that the nor- 
malized nucleation density 

M(/,e) = n(/,e)/n(/,l) (8) 

obeys an unusual logarithmic scaling law of the form 

lnA'f(/,e) ~ (In/) .g(lne/ln/). (9) 

This scaling form was also proposed by Kadanoff et al. 
and Tang in the context of multiscaling in self-organized 
criticality j|, |l^. As shown in Fig. |^this scaling form 
leads to a convincing data collapse. More recently the 
same type of scaling has also been observed in one- 
dimensional systems, so that we can rule out logarith- 
mic corrections at the marginal dimension of random 
walks as the origin for this type of non-conventional seal- 
ing (l|,|l|. 

In Fig. ^the left terminal points of the curves at 17(1) = 

0. 44(2) can be used to determine 7 as follows. It is known 
that the nucleation density after the first pulse grows 
linearly with the intensity jSj while the nucleation density 
after completion of one monolayer grows as n(/, 1) oc 
I"^'^ (see Eq. (^). Therefore, the normalized nucleation 
density after the first pulse scales as Af (/, /) oc /i-^i^ 
that 

1-47 

g[l) = l-2v=^—^. (10) 

Solving this equation we obtain 7 = 0.179(4) in agree- 
ment with the estimate in Fig. ^. 

In order to understand how this unusual scaling be- 
havior for PLD crosses over to the ordinary power-law 
scaling of MBE let us now turn to the case of finite D/F^ 
where the system is characterized by a typical length 
scale €0 ''^ {D/FY/^. Generalizing the results of Ref. |^ 
we consider islands with an arbitrary fractal dimension, 

1. e., we regard 7 as a free parameter. 
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FIG. 5: Data collapse according to the scaling form 
The inset shows a double-logarithmic plot of In A// In/ vs. 
In O/ In 7. The line is slightly curved with local slopes vary- 
ing from 2.4 to 2.0. 



In MBE the nucleation density is known to exhibit 
ordinary power law scaling of the form 



(11) 



where / is a scaling function with the asymptotic behav- 
ior O 



for < z < 1 
2:1/3 for 1 < z < 



(12) 



The upper bound z^ax is determined by the condition 
that the whole surface is covered by islands so that no 
further nucleation in the respective layer is possible 0. 
Because of Ic ~ {D / FY^~^ the length scale is related 
to the crossover intensity by 



^ 7-1/(4-87) 

Using this scaling form we obtain the expression 

3-1/3 



for < z <C 1 

for 1 < Z < Zniax , 



(13) 



(14) 



where z — QIc ^^^'^ ^'^^ ■ 
trapolating to z — i\ we arrive at 



Taking the logarithm and ex- 



In M(Je,e) 
In/. 
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^ In /e 

Ine 
3 In 7c 



67-3 



for 



1 ^ Ine 
2-47 ^ In/c 

In e 



for « Tin; « 2-47 



(15) 



In the limit / ~ ^ (i.e., D/F ^ 00) the crossover 
between both regimes becomes sharper and converges to 
a piecewise linear curve, which is shown in Fig. ^ as a 
dashed line. The crossover point is located at 



InGe/ln/e 
lnAf(/„e)/ln/e 



1/(2-47), 
1/(6-127), 



(16) 



Surprisingly, the crossover point lies on the collapsed 
curves for PLD within numerical errors. This is plau- 
sible for the following reasons. On the one hand the 
PLD-curve must be an upper bound for the MBE-curve 
since the island density in PLD is always larger. On 
the other hand, if the gap between the two extrapolated 
curves did not close at the crossover between PLD- and 
MBE-behavior, it would imply that there is an additional 
characteristic length in the system, for which we have no 
evidence. 

Furthermore, we note that the scaling function itself 
roughly follows a power law 



g{z) ~ g(l)z 



/3 



(17) 



As shown in the double-logarithmic inset of Fig. ||, the 
effective exponent 13 varies between 2.4 and 2.0. In or- 
der to verify this estimate we derive the exponent (3 from 
Eq. (|l^) by assuming that the crossover point lies ex- 
actly on the collapsed curves for PLD. This leads to an 
expression for (3 in terms of the exponent 7, namely 



ln(6 - 247) 
ln(2 - 47) 



2.4 for compact islands 

2.15 for DLA fractals ' ^ ' 



Since this is just the range in which the numerical values 
for (3 vary, we are led to the conclusion that the effec- 
tive fractal dimensions of the islands for low and high 
coverages are different. For low coverages the islands are 
spaced relatively far apart so that the growth is DLA-like 
while for high coverages the islands coalesce and become 
more and more compact. 



IV. INFLUENCE OF EHRLICH-SCHWOEBEL 
BARRIERS 



In most experimental situations interlayer transport is 
reduced by an additional energy barrier i?ES which the 
atoms have to overcome when hopping down an edge of 
an island |l4| . This barrier is called Ehrlich-Schwoebel 
(ESVbarrier and is of the order of about 0.1 eV for met- 
als . A useful measure for the influence of this barrier 
is the Schwoebel length £es , which is defined as n& 



£es = a- exp 



kBT 



(19) 



where a — 1 denotes the lattice constant. For MBE 
it is known that Ehrlich-Schwoebel barriers impede in- 
terlayer transport. This leads to a growth instability 
which was predicted and first theoretically investigated 
by Villain |l^ followed by others ||l^ and has also been 
observed experimentally ||l^. In contrast to MBE, not 
much is known about the influence of Ehrlich-Schwoebel 
barriers on PLD. However, it has been observed exper- 
imentally that PLD leads to better growth results than 
MBE in certain situations where Ehrlich-Schwoebel bar- 
riers are present 0. 



The influence of ES barriers on PLD is twofold. On 
the one hand, for high intensities many adatoms are de- 
posited on the same island, which should increase the 
influence of ES barriers. On the other hand, in PLD the 
islands arc much smaller so that adatoms tend to leave an 
island very quickly, thereby reducing the influence of ES 
barriers. Therefore, the question arises whether the ex- 
perimentally observed improved quality of layer-by-layer 
growth in PLD can be related to a reduced influence of 
ES barriers and whether it is possible to choose the pa- 
rameters of the model in such a way that PLD produces 
smoother surfaces than MBE. 

This question has been addressed previously by 
Schinzer et al. |Q. They studied the special case where 
one intense pulse (0.23 ML) is deposited at the begin- 
ning of each monolayer, while the remaining atoms are 
deposited with a continuous flux. The aim was to check 
an idea of Rosenfeld et al. |2l| , that it should improve the 
smoothness of the surface, if the islands are forced to be 
much smaller than the diffusion length of the atoms de- 
posited after the pulse. Surprisingly, Schinzer et al. pc| ] 
found that this is not the case in their simulation. On the 
contrary, in the presence of the pulse the surface became 
rougher. 

A key concept for describing the influence of ES- 
barriers is the time it takes before the first nucleation 
in the second layer takes place. Layer-by-layer growth 
requires that nucleations in the second layer do not start 
significantly before the first layer is completed jl^. In 
the case of MBE the study of the second-layer nucleation 
time turned out to be very useful in order to predict 
the growth mode for a given set of parameters p3[ . 
In the following we demonstrate that this concept can 
be successfully applied to PLD as well. To this end 
we first investigate the growth behavior of PLD in the 
limit of D/F — > oo. It will be shown that for any finite 
Ehrlich-Schwoebel barrier MBE produces a better layer- 
by-layer growth than PLD in agreement with the finding 
of Schinzer et al.[|o). Then we compare MBE and PLD 
for different values of / and D/F over a wide range of 
Ehrlich-Schwoebel barriers. It turns out that there is in- 
deed a regime where PLD produces better growth results 
than MBE. 



0.15 




20 30 
Island Area 



0.1 



0.05 



1 ' 

(b) 


1 1 1 




1=0.01 M"'- 








1=0.005" 

1,1 



10 



20 30 
Island Area 



40 



50 



FIG. 6: Numerical measurement of the second-layer nu- 
cleation probability compared to the theoretical prediction. 
After each pulse pnuc is estimated as the fraction of islands 
with area A, where a second-layer nucleation happened for 
the first time, averaging over several runs. The island area 
is restricted to A < 50, where coalescence does not yet play 
a role, (a) The second layer nucleation probability pnuc for 
I = 0.01 versus island size A for different Ehrlich-Schwoebel 
barriers in the limit D/F oo. From bottom to top the 
Schwoebel lengths are £es = lO", 10\ 10^ 10^ 10*, 10^ (b) 
Plot of Pnuc versus A for two different intensities and for high 
Ehrlich-Schwoebel barriers ^es = 10^ lO'*, 10^ 



A. Infinite D/F 



We start by explaining why in the limit D/F — > oo 
PLD does not delay the Villain instability in spite of the 
island size reduction compared to MBE, provided that 
the ES-barrier is sufficiently high and not leaky e.g. at 
kink sites. In this limit the time scales are separated as 
follows. Because of the high ES barrier the nucleation 
time of adatoms on top of islands is much smaller than 
their residence time, while for D/F oo the residence 
time is in turn much smaller than the time interval be- 
tween two pulses. This means that we can restrict our 



analysis to a single pulse. 

The second layer nucleation probability pencil , A) is 
defined as the probability that there is at least one nu- 
cleation event on islands with area A after deposition of 
a pulse with intensity /. For high ES barrier this means 
that during a single pulse at least two atoms have to be 
deposited on the same island. Obviously, the probability 
of depositing k atoms on an island with area A during 
one pulse is given by a Poisson distribution with average 
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FIG. 7: (a) Squared surface width versus time for MBE and 
PLD for different intensities. Tlie growth parameters are 
D/F — 10*^ and ^bs = 1- (b) Squared surface width versus 
intensity for D/F = 10* and ^es = lO"*. 



Pk 



k\ 



(20) 



Hence the nucleation probability can be expressed as 



P„uc(/,^) = = l-e"'^ (1 + M) 

~ I^A^ + 0{I^A^). (21) 

Fig. ||a shows the numerical results for pnuc for / — 0.01 
and different ES-barriers. As can be seen, for increasing 
ES barrier the measured curves approach the predicted 
one, (|l]), which is shown as a dashed line. The agree- 
ment is good for ^ 10'^, where the nucleation time 
is much shorter than the residence time. However, the 
agreement is not convincing for island areas A < 10, 
where the discrete lattice spacing starts to play a role. 
Moreover, the upper curves deviate for A > AO due to co- 
alescence of large islands. This explanation is supported 
by the results shown in Fig. where two different in- 
tensities are compared. Indeed, for the lower intensity 



/ = 0.005, where coalescence starts at larger island sizes, 
the measured and the predicted curves agree much bet- 
ter. Thus we conclude that the second layer nucleation 
probability pnuc is adequately described by (|2l|). 

With Eq. (|2^) we can now answer which of the oppos- 
ing trends in PLD - island size reduction and increasing 
nucleation probability on top of islands - will dominate. 
In the limit D/F —>■ oo, where the critical intensity Ic 
tends to zero, the island area scales as A ^ j-2i^ 
for all intensities. Together with Eq. ( ^ ) one obtains to 
leading order 



Pn 



oc 



/2(1- 



2u) 



(22) 



Using the previous estimate v 2± 0.26 the exponent is 
given by 2(1 — 2^) ~ 1. Therefore, the second layer nu- 
cleation probability Pnuc grows with increasing intensity, 
enhancing the Villain instability. Contrarily, MBE with 
infinite D/F always shows perfect layer- by-layer growth, 
even for high but finite ES-barriers. Thus, we conclude 
that in the limit D/F — > oo PLD cannot improve layer- 
by-layer growth. In order to apply this result to practical 
situations, one has to find the lower boundary for D/F, 
below which the behavior is different. This question will 
be addressed now. 



B. Comparison of PLD and MBE for finite D/F 

For perfect layer-by-layer growth the squared surface 
width w'^ = {{h/ — {h'^)) is known to oscillate between 
zero for completed monolayers and 1/4 for half mono- 
layers . Without ES barriers such oscillations can be 
seen in MBE as well as in PLD. As shown in Fig. they 
are most pronounced in MBE while in PLD they become 
more and more damped as the intensity increases. More- 
over, it can be seen that in PLD without ES barriers the 
roughness always increases with increasing pulse inten- 
sity 

Fig. 1^ shows the corresponding result for a very high 
Ehrlich-Schwoebel length ^es = 10^. There are no os- 
cillations since the surface roughens very quickly due to 
the Villain instability . Surprisingly, the width is now 
maximal for MBE and decreases with increasing pulse in- 
tensity, i.e., the trend is reversed. This reversal can even 
be observed visually by monitoring the interface at dif- 
ferent intensities, as shown in Fig ^ 

Let us now investigate this reversal in more detail. 
Since the oscillations in PLD are extremely weak or not 
present at all, a numerical measurement of the damp- 
ing time (such as in Ref. for MBE) is not feasible. 
Instead we monitored the squared surface width after de- 
position of a two monolayers - a coverage which is also 
relevant for experimental applications Although the 
choice of two monolayers is arbitrary and does not per- 
mit a rigorous quantitative analysis, this criterion is very 
simple and tells us for which parameters PLD produces 
smoother surfaces than MBE. 



7 




FIG. 8: 1+1-dimensional PLD after deposition of 20 ML for different pulse intensities. The nucleation sites are shown as 
white dots. Left: Without Ehrlich-Schwoebel barrier the roughness increases with increasing intensity. Right: If the Ehrlich- 
Schwoebel barrier is high enough, the autocorrelation among nucleation sites is more pronounced than for ^es = 1, and the 
roughness decreases with increasing intensity. 



The roughness after deposition of two monolayers is 
plotted in Fig. ^. In agreement with previous results 
for MBE all curves increase monotonously, i.e., the 
roughness increases with increasing Ehrlich-Schwoebel 
barrier. However, the curves for MBE and PLD cross 
each other. The typical ES barrier, where this crossing 
takes place, varies roughly as £es ~ (-D/F)^/^. This ob- 
servation is in agreement with our previous result that for 
D /F oo surfaces grown by MBE are always smoother 
compared to PLD. 

Thus there is indeed a range of parameters where 
PLD with a high intensity produces smoother surfaces 
than MBE, even if the atoms are deposited with ther- 
mal energy. However, the Ehrlich-Schwoebel barriers 
are unphysically large (typical experimental values are 
£es ~ 1...10). Therefore we believe that nonthermal 
energy deposition effects are important for explaining the 
experimentally observed growth improvement in PLD. 

Nevertheless it is interesting to address the question, 
what causes the crossing of the curves in Fig. ||. As we 
have explained, second layer nucleation becomes more 
likely the larger the pulse intensity is in spite of the de- 
creasing island size. Therefore one would always expect 
the roughness to increase with pulse intensity, as it is 
indeed the case for sufficiently small Ehrlich-Schwoebel 
barriers. For €es ^ however, there is no interlayer 
transport and the roughness is the same as in random 
deposition, irrespective of the pulse intensity. Thus, the 
curves in Fig. |^ are expected to saturate eventually at 
the same roughness. Thus the observed reversal results 
from an interplay of a large but finite ES barrier and a 
finite value of D/F. 

To understand the reversal let us consider how the Vil- 
lain instability unfolds for large but finite ES barriers. 
Depending on the intensity, the first pulse leads to the 
formation of many small islands. After the first pulse 
there are two temporal regimes of different roughening 
behavior. Because of the Villain instability effective up- 




FIG. 9: The squared surface width at B = 2ML versus the 
Schwoebel length for MBE and PLD with various intensities 
and for D/F = 10*. The curve for MBE is dashed. 



hill currents lead to a quick formation of large mounds. 
As new nucleations preferentially take place on top of 
these mounds, the nucleations are vertically aligned , 
as shown on the r.h.s. of Fig. [VB. Therefore, the num- 



ber and the lateral size of the mounds is essentially deter- 
mined by the initial configuration of small islands after 
the first pulse. This determines at least for transient 
times also the height of the mounds and therefore the 
roughness of the surface. In our simulation model we 
do not observe slope selection. In a more realistic case 
one would expect that the mounds grow until their edges 
reach a critical slope where the uphill current becomes 
zero. Then the process enters a different temporal regime 
where the mounds compete with one another, leading to 
an effective coarsening process where the roughness in- 
creases only slowly. 

By increasing the pulse intensity second layer nucle- 
ations become more likely so that the Villain instabil- 



ity is accelerated. However, the typical roughness from 
where on the process enters the second temporal regime 
of slow coarsening depends mainly on the initial density 
of islands after the first pulse. Therefore, this typical 
roughness decreases with increasing intensity. After suf- 
ficiently long time this effect dominates and leads to the 
observed crossings of the curves in Fig. ^. Obviously, this 
mechanism works in any dimension. 



V. CONCLUSION 

In the present paper we have studied a model for pulsed 
laser deposition with and without Ehrlich-Schwoebel bar- 
riers. The model assumes that the particles are de- 
posited at thermal energies where the transient mobility 
of adatoms as well as implantation effects can be ne- 
glected. First we investigated the case without Ehrlich- 
Schwoebel barriers. For low pulse intensities PLD dis- 
plays essentially the same behavior as MBE. Increasing 
the intensity it crosses over to a different behavior charac- 
terized by the nucleation of many small islands after the 
first pulse. The behavior in this regime can be analyzed 
by studying the time-dependent nucleation density. Ex- 
tending previous results we have shown that in the limit 
D I F — > oo this quantity displays an unusual type of log- 
arithmic scaling behavior. 

Turning to PLD with Ehrlich-Schwoebel barriers with- 
out leakages we first showed that in the limit D/F ^ oo 
the second layer nucleation probability grows with in- 
creasing intensity, enhancing the Villain instability. This 
means that in this limit PLD cannot improve layer-by- 
layer growth. For finite D/F, however, the situation is 
different. Studying the surface width for the Schwoebel 
lengths between 1 and 10"* we found that PLD produces 
a smoother surface than MBE if the barrier is strong 
enough. This reversal can be explained by the influence 
of the initial nucleation density after the first pulse on the 
roughening due to the Villain instability. However, we 
believe that this mechanism alone cannot explain the ex- 
perimentally observed improvement of layerwise growth 
in PLD where the barriers are much lower. In order to 
describe this effect more accurately the transient mobility 
of adatoms would have to be taken into account. 

As we pointed out, Schinzer et al. |po| had reached 
a similar conclusion that pulsed deposition leads to in- 
creased roughness. However, this statement, which con- 
tradicts an idea by Rosenfeld et al. cannot be uni- 
versally valid, as our more detailed investigation showed. 
In our opinion the negative effect of pulses is mainly due 
to the fact that the probability of multiple deposition on 
top of islands is increased. This indicates, that the very 
high intensity of the single pulse per monolayer was re- 
sponsible for the roughness in the simulations of Schinzer 
et al. @. 

However, Schinzer et al. also discovered a way to 
improve MBE with pulsed deposition: If desorption can- 




D/F 



FIG. 10: Island distance for MBE and PLD with J-L^ = 1 as 
a function of D/F for a 400 x 400 lattice. The island distance 
is measured at 0.2ML. 

not be neglected, atoms will preferentially evaporate from 
the top terraces. Thereby the negative effect of multiple 
deposition can be compensated, so that one gains the 
benefit of enforcing small islands by the pulse in the be- 
ginning of the monolayer. 

Our model describes only a limiting case of pulsed 
laser deposition. It has the virtue of making new scaling 
concepts clear. In practice the atoms are seldomly de- 
posited with thermal energy in pulsed laser deposition. 
In particular it would be interesting to work out the 
effect of transient mobility due to heating by the pulse. 
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APPENDIX A: FINITE-SIZE EFFECTS IN THE 
LIMIT OF LOW INTENSITY 

As discussed in Sec. II, for small intensities PLD dis- 
plays essentially the same growth behavior as MBE. How- 
ever, on finite lattices both processes are not exactly 
equivalent. This can be seen in Fig. |l^, where the av- 
erage island distances for MBE and for PLD with an 
intensity of one particle per pulse / = are compared 
over a wide range of D/F. For very low values oi D/F 
the curves coincide but do not follow a power law because 
of lattice effects. Only in the range 10^ < D/F < lO'^ 
the curves follow approximately the expected power law 
Id oc {D/F)''. Finally, for D/F > lO'^ the island dis- 
tances of MBE and PLD differ increasingly from each 
other. 

In order to demonstrate that this discrepancy between 
PLD and MBE is a finite size effect we determine the 
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FIG. 11: Values of D/F where the island distances in PLD 
with IL^ = 1 and in MBE differ at least by a certain per- 
centage specified in the legend. The slopes of the power-law 
fits from bottom to top are estimated hy x = 3.0(2), 3.1(2), 
3.1(2), and 3.2(2). 

values of {D/F)dov from where on the island distances 
in both models differ by a certain factor. As shown in 
Fig. |ll| these values increase algebraically with the sys- 
tem size as 

(i?/F)dev K , (Al) 



where the exponent x is found to be close to 3. This 
power law behavior can be explained as follows. In both 
cases single atoms are deposited and the average time 
between two depositions t = 1/L^ is the same. However, 
in PLD the deposition takes place in constant time inter- 
vals, whereas in MBE atoms are randomly deposited so 
that the time intervals between deposition events obey 
a Poissonian distribution P(t) = L^e~^ with the vari- 
ance a = \J (r^) — (r)2 = 1/L^. If such a fluctuation 
leads to a r which is smaller than the average time in- 
terval between two depositions, more nucleations will be 
produced. As the formation of nucleations is irreversible 
this enhancement is not compensated by fluctuations in 
opposite direction where r is large. Therefore, the influ- 
ence of fluctuations increases the number of nucleations, 
leading to a smaller average island distance than in PLD. 
However, this effect can only be seen if the fluctuations 
are strong enough, i.e., they have to be at least of the 
same order as / {D / F), which is the average diffu- 
sion time (in ML) before an adatom reaches the edge of 
an island or another adatom. Since to ^ (D/F)'* this 
argument leads to the scaling relation 

(i?/i^)dev cx L\ X=T^- (A2) 



For 7=1/6 the exponent is x = 3 while for DLA it is 
given by x — 3.1. This result is in fair agreement with 
the numerically determined exponents in Fig. 
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